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Wave Motion in Hydrodynamics. 

Br A. G-. GtKEenhill, Woolwich, England. 



The mathematical subject of Hydrodynamics is still, in some respects, in an 
elementary stage, insomuch as every fresh problem solved constitutes a distinct 
advance of the subject. 

One of the most important applications of the theory of Hydrodynamics is 
to the question of the motion of Waves under gravity and other causes, and as 
the investigations on this subject are for the most part scattered about in various 
scientific periodicals, I propose in this article to collect together the chief results 
hitherto obtained, and to give also a general connected account of the mathe- 
matical theory, at the same time attempting to develop it in some directions. 

In the mathematical treatment of Wave Motion we are constrained at 
present to employ the approximation of supposing the velocities of the liquid 
particles due to the wave motion to be sufficiently small for the squares, etc., of 
the particle velocities to be neglected ; although it is singular that this approxi- 
mation is not required in the first problem of wave motion ever solved, discov- 
ered by G-erstner in 1802, and afterwards independently by Rankine in 1862 
(Stokes, Mathematical and Physical Papers, I, p. 219). 

A list of the principal papers on the subject of Wave Motion and of their 
authors will be found in the Report on Recent Progress in Hydrodynamics, 
by W. M. Hicks, F. R. S., presented to the British Association. 

1. The most convenient order to employ in the mathematical treatment of 
a problem in the subject of Wave Motion is: (I) The determination of the 
velocity function <|> , or stream function tJ/, satisfying the equation of continuity ; 
(II) The determination of the boundary conditions to be satisfied at the sides of 
the containing vessel ; (III) The most difficult part, the determination of the 
conditions to be satisfied in order that the free surface should be a surface of 
equal pressure, or, more generally, at the surface of separation of two liquids 
there should be no discontinuity of pressure. 
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To secure uniformity of notation in the treatment of waves under gravity, 
we shall suppose the co-ordinate axis Oz drawn vertically upwards, and the 
plane xOy taken generally in the undisturbed horizontal plane of the surface of 
separation of two liquids at which the wave motion is apparent, and then the 
axis Ox will be taken in the direction of propagation of the waves when straight- 
crested, and, therefore, perpendicular to the crests. 

Then, to determine the wave motion in still water (still except for the 
slight disturbance of the wave motion), we must determine a velocity function 
$ . (I) Satisfying the equation of continuity for liquids : 

av .**.?*_ 

dx 1 ^ dtf ^ dz 2 
(II) Satisfying the condition 

-f£ = 

av 

at a fixed boundary, dv denoting an element of the outward-drawn normal of 
the boundary, or, more generally, 

— — L ^Z f) 

dv n ' 

the normal component velocity of the boundary, when movable. (Ill) Satisfy- 
ing at the surface 2=0, supposed a surface of equal or of no discontinuity of 
pressure, the dynamical equation 

* + „+£=* 

a constant, neglecting the squares of the velocities of the liquid particles. 

At a free surface p is constant, and, therefore, -£- = , so that, 37 denoting 
the elevation of the free surface, 

9 dt + dt 

Bufc W = ~dz ' 

and ^ = -J^<p, 

where I denotes the length of the equivalent simple pendulum of the wave 

motion, so that 7 da> 

l ~dz--*> 
at the free surface 3=0, or z= h, some constant. 
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2. When surface tension T is taken into account, this equation must be 

replaced by 7^4. T[_ <?$_ — A 

dz go dz 6 

an equation due to Kolacek {vide Fortschritte der Mathematik, 1878). For, if 
dp denotes the excess of pressure in the liquid just below the capillary film over 
the external pressure above, 

& + ^ + _|£ = o. 

p at 

But r and r 1 , denoting the radii of curvature of any two vertical sections of the 
free surfaces by perpendicular planes, 

to one order of approximation, so that 

and differentiating with respect to t, and replacing -p- by -SjE- , remembering 

dt dz 

also that c^> , 3> _ _ 3> 

dx> + dy % dz' ' 

3. Digression on the Hyperbolic Functions. 

In the course of our investigations we shall require certain functions, called 
Hyperbolic Functions, from their connection with the hyperbola, which are 
analogous to the functions of the circle defined in ordinary trigonometry. As 
these functions are not defined and explained in all the ordinary text-books, we 
shall, for convenience, proceed to do so as follows : 

(I) — (e* + e~*) is called the hyperbolic cosine oiv, and is denoted by cosh v. 

(II) -o~(e* — e~ v ) is called the hyperbolic sine of v, and is denoted by sinh v. 

(III) „_. _ 8 = — r- is called the hyperbolic tangent of v, and is denoted by 
tanh v ; and so on, by analogy, with the rest of the circular functions. 
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From the exponential values of the cosine and sine, viz., 

cos u = — (e iu + e~ tw ) , sin u = — r (e iw — e~ iu ) , 

when i denotes V — 1, we see, by putting u = iv, that cos iv = coshv, 
sin iv = i sinh «, tan iv = i tanh », etc. ; also, 

cos (u + »») = cos u cosh « — i sin w sinh v . 
sin (u + *t>) = sin u cosh « + i cos m sinh v, etc., 
formulae of great use hereafter. Therefore, also, 

cosh (u + iv) = cosh u cos t? — i sinh u sin « , 
sinh (u + iv) = sinh u cos » — i cosh w sin « . 

Analogous to the ordinary formulae of circular trigonometry, we have 

cosh (a + /?) = cosh a cosh ft + sinh a sinh ft , 
sinh (a + ft) = sinh a cosh ft + cosh a sinh ft, 

, , . /1N tanh a + tanh fi . 
tanh < a + ^ = i + tanh « tanh /9 ' 

sinh y + sinh h = 2 sinh — (7 + 5) cosh — (7 — 8) , 



sinh 7 — sinh 5=2 cosh — (7 + 5) sinh -„- (7 — 5) , 
cosh 7 + cosh 5 = 2 cosh -— (7 + 5) cosh -5- (7 — 5) , 
cosh 7 — cosh 5=2 sinh — (7 + 5) sinh — (7 — 5) , 



and so on. 



4. TPaves in Still Water of Uniform Depth. 



Supposing straight-crested waves of length X propagated in the direction 
of the axis of x with velocity U, we may begin by supposing the velocity 
function <p = f(z) cos (mx — nt) , 

where w = 27t/A, n= 2x17/1, and n % = g/l, I denoting the length of the 
equivalent simple pendulum. 

Then, from the equation of continuity, 

f— y=o, 

the solution of which is 

f(z) = ae m ' + ie- ms ; 

VOL. IX. 
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or, using hyperbolic functions, 

f{z) = P cosh mz + Q sinh mz ; 
or, subject to the condition that 

4£=0, whenz=— h, 

oz 

h denoting the depth of the water, 

f(z) = A cosh m (z + h) ', 
so that A = A cosh w (z + h) cos (mcc — nt) . 

Then, at the free surface z = , 

or, m? sinh wi/i = cosh wi/i ; 

or, <ml= cothw/t. 

Then, 7J2 = 4=^y 

= -2- tanh w/ = -2— tanh— — , 

the well-known expression for the wave velocity. 

When h/X is small, we can replace tanh (2nh/X) by 2tHi/%, and then 

U* = gh, 
Kelland, Scott, Eussell and Green's expressions for the wave velocity when the 
wave length is great compared with the depth of water. 

When h/X is large, we can replace tanh (27th/ X) by unity, and then 

U* = g%/27t, 
agreeing with Gerstner's and Rankine's expressions for the wave velocity in 
water of great depth. 

5. Next, suppose there is a surface tension T at the free surface ; then the 

condition 7 9<2> , Tl d% , 

dz go dz* 
when z = , leads to the relation 

Tl 
ml sinh mh A m 3 sinh mh = cosh mli ; 

99 



TV— n 

or, U" = — » 



q( 1 } tanh ?n7i 

Ho> gp J 



99 
pX 



/ gX , 2rtT\ , , 2/rh 
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the general expression for the wave velocity under gravity and surface tension 
combined. 

When 2, is large or T small, we may put 

TT „ gk 1 , 2nh 

as above ; but when X is small, the term g^/2n is insensible, and we can put 

TT „ 2xT, 27th 
U % = — r tanh — - , 

pA A 

the velocity of propagation of ripples of wave length % due to surface tension T. 
Supposing the depth of water h sufficiently large for tanh (27th/ X) to be 
replaced by unity, then 

m =■§; + jr = 2 W + w i ~ \/~pn ' 

so that the minimum value of £7 is 

and the X = 2%s/ (T/go). 

Sir W. Thomson proposes to distinguish by the name of ripples those 
waves whose length is less than the above critical value of a. ( Phil. Mag. (4) 
xlii). 

6. A slight extension of this problem may be made by supposing the 
capillary film of the surface to be replaced by a flexible cloth of uniform tension 
^and uniform superficial density a resting on the surface of the liquid. 

Then, assuming for the liquid motion a velocity function 

q> = A cosh m (z + h) cos (mx — nt) , 

as before, supposing A a small constant factor, and denoting by ri the elevation 
of the surface, then, when z = , 

~ = -J±- = m A sinh mh cos (mx — nt) , 
at dz 

and, therefore, x> = -A sinh mh sin (mx — nt). 

n v 

Denoting by dp the excess of pressure just below the cloth over the atmos- 
pheric pressure above, then 

3? j_ „„ j. 3<£ 
P 

or, dp — go — n A sinh mh sin (mx — nt) + n A cosh mh sin (mx — nt) = . 



+ W + 3 = 0i 
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But the equation of motion of the cloth is 

so that, dropping the common factor A sin (mx — nt) , 

cmn sinh mh = T — sinh nth + go — - sinh mh — on cosh mh; 

or U 2 = — = ^ /m + Tm 

m a /? coth mh + <"** ' 

giving U the velocity of propagation of the waves, and reducing, when a= , 
to the preceding case of a capillary film. 

7. Waves in Ice of Uniform Thickness Besting on Water of Uniform Depth. 

If the water is covered with ice, then the equation of vibration of the 
surface must be replaced by 

where L denotes the flexural rigidity of the ice, the vibrations being now of the 
nature called lateral vibrations (Rayleigh, Theory of Sound, I, Chapter 8, 
§163), the inertia of each vertical section of the ice being supposed concentrated 
at the centre. 

Then if e denotes the thickness of the ice, and E Young's modulus of 
elasticity, r 1 ,„ n 

J 12" ' a ~ e ?' 

supposing the ice of the same density as water, so that now 

amn sinh mh = — L — sinh mh + go — sinh mh — on cosh mh ; 
or, U % = — = 9P/ m + Lm * 



p coth mh + am 

g%/2n+ y^^E/o% 3 
~ coth(2^A/>l) + 2w6/A ' 
giving the velocity of propagation of waves of length X in ice of thickness e, 
resting on water of uniform depth h . 

It is remarkable that ice was the first substance for which an experimental 
determination of E was attempted, as described in Young's Lectures on Natural 
Philosophy. 
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8. Waves in Water of Uniform Depth Established and Maintained by Impinging 

Waves of Sound in Air. 

Suppose the preceding kind of wave motion in conjunction with plane 
waves of sound impinging at an angle (3, we have thus an illustration of a 
forced, or rather controlled, wave motion in the water due to free waves in the 
air. Let £ = B sin { m (x sin /? — acos/3) — nt-\-a) 

represent the normal displacement in the incident wave of sound, and 

£i = B x sin { m (x sin /? + z cos (3) — nt + a x \ , 
in the reflected wave ; and let 

vi = b sin {mx sin /? — nt) 

represent the displacement of the surface of the water. Then, at this surface, 
we must have n = (& — £) cos (3 , 

when 2=0, in order that there should be no separation of the air from the 
water, so that 

b sec/? sin (mx sin l3~nt) = B 1 sin (ma: sin /3— nt + aj — B sin (mx sin/3 — nt-\- a) 
for all values of x, leading to the equations 

/?! cos (*! — B cos a = b sec (3) 
/?! sin a x — B sin a = 
The velocity function of the motion in the water must be of the form 
<p = A cosh [m(z + h) sin/?}- cos (mx sin (3 — nt); 



\- 

0) 




/////////////////////////////////////////////////////////////////////////////// 



and then, since 

when z = , therefore, 



— n& = mJl sin (3 sinh (w7i. sin (3) . 
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From the dynamical equation 

f +«■+!?=* 

we have, denoting by dp the periodic part of p just below the surface of the 
Water ' f +a « + f=0; 

or, dp = — goYi — o-^- 

= — \gpb + on A cosh (mh sin /3) } sin (ma; sin (3 — nt) . 

Suppose, now, at the surface of separation o*f the air and the water there is 
a film or cloth of tension T and superficial density a ; and we might also, if we 
like, suppose the film to possess flexural rigidity L like ice, without much addi- 
tional complication ; then, at the surface of separation, 

when z = , where dp' denotes the periodic part of the pressure in the air due 
to the wave motion. 

Now, o' denoting the density of air and a the velocity of the sound waves, 
we have ri* t p' 

and the cubical elasticity 



m s 'pt ' 




/ d P' i 
,dp' , 

yp ^r= —yps — 

r 




— a\'i 



so that dp ■ 

s denoting the cubical expansion ; so that 

s = Tx ( * 1+ * /sin ® + 1 & ~~ ® cos ^ 

= mB x cos { m (x sin (3 + z cos (3) — nt + a x \ 
+ mB cos \m(x sin (3 — z cos (3) — nt -\- a\. 
Therefore, when a = 0, 

dp' = — cPo'rn \ B x cos (nix sin (3 — nt -\- ai) -\- B cos (mx sin /3 — «£ + a) \ 
= a s p'm (B x sin aj -f- -B sin a) sin (mx sin /3 — rc£) , 
provided that B 1 cos a x + B cos a = , 

which, combined with the previous equations, gives 

Bi sin a _ cos a 

B sin a! cos a t ' 
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or, tan a = — tan a x , a = — a x , 

and B 1 = — B = •— 5 sec a sec /? . 

Then, since 6 = — — A sin {3 sinh (wi# sin /?) , 

the boundary condition becomes, when the common factor A sin (mx sin (3 — nt) 
is dropped, 

cwm sin {3 sinh (mjA sin (3) = T — sin 3 (3 sinh (mh sin /?) 



— L — sin 5 ft sinh (mh sin /?) + gp — sin (3 sinh (mh sin /?) 

— pn cosh («iA sin /3) — a % p' — tan a tan (3 sinh (»i& sin (3) ; 
or, p'n tan a tan /? = — omn sin /? + T 7 — sin 3 /? 

— L — sin 5 /? + gp — sin /? — pn coth (mA sin /?) , 

giving tan a , and therefore a , the change of phase of the sound wave in being 
reflected at the surface of the water. 

9. Reflection and Refraction of Plane Waves of Sound by a Plane Curtain. 

Suppose, now, that the cloth, instead of resting on the surface of a liquid, is 
the plane surface of separation of two elastic fluids of different densities p and p', 
but necessarily of the same pressure p when at rest, the cloth being now 
supposed vertical to abstract the curving effect of gravity; let us now inves- 
tigate the reflection and refraction of plane waves of sound, impinging and 
being reflected at an angle 8 in the first medium of density p , and being refracted 
at an angle (3' in the second medium of density p'. 

Denoting as before by 

£ = B sin \m (x sin (3 — z cos (3) — nt + a \ 
the normal displacement of the incident waves, and by 

& = B x sin \m (x sin (3 + z cos 8) — nt + ocj J- 
of the reflected, and by 

% = B' sin \m'(x sin (3' — z cos 8') — nt + a' } 
of the refracted wave ; then 
(I) m sin (3 = m! sin 8', the Law of Refraction. 



72 



G-KEENHILL : Wave Motion in Hydrodynamics. 



Also, denoting the displacement of the surface of separation, the cloth or 
curtain, by rf=b sin {nix sin /? — nt) ; 

then, as before, when 2 = 0, 

V = (£1 — £) cos /3 = — £' cos /3', 



leading to 




B x cos a x — B cos a = sec /3 , 
J?! sin a x — B sin a = , 

B' cos a' = — b sec /3' , 

B sin a' = ; 

so that a! = , and 5' = — b sec /?'. 

For the motion of the curtain, 



cfty 



S 



where 3p = — a?ps, dp' =■ — a' 2 p's', 

s and s' denoting the cubical expansion in the two media in the neighborhood of 

a = ; also, a and a' denoting the velocity of sound in the two media. Then 

s = mB x cos (mx sin (3 — nt + ai) + inB cos (ma; — «£ + a) , 

s' = m'B cos (w2jc sin /? — nt) ; 
so that 

— c -p- + Sr-rl = (<7W 2 — Trtf sin 2 /3) 6 sin (mx sin /3 — nt) 

— dp — dp' = — a 2 pm \ B x cos (mx sin /? — n£ ■+- c^) + B cos (mx — nt + a) } 
+ a m p'm'B' cos (mx sin /? — «£) , 
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for all values of x and t, leading to the conditions 

(an 9 — Trn % sin 2 ft) b = a?om (B t sin a^ + B sin a) 

and = — a % om (B x cos a t -\- B cos a) -f a' 2 pW.B'. 

Therefore, since a 2 = « 2 /w 2 , a' 2 = rc 2 /m' 2 , 

„ . „ . <m 2 — Tm 2 sin 2 i3 7 <ra 2 — Tm 2 sin 2 £ T 

if, sin ai = -d sin a = ^ — 5 L mo = —5 b , 

2pn 2 2d i pm 

and B, cos a, + -B cos a = „' i?', 

aym 

= -^6 sec/3'; 

also, J5 X cos a! — -B cos a = J sec /? ; 

, , , „ crW sec fl — a'Vm' sec i? 7 
so that B x cos a x = — —-^ — - b , 

„ d l pm sec /3 + a'^p'm' sec # T 

5cosa= - 4r ! 2 — ■ b , 

2a i pm 

■ . , a % pm sec 8 4- a' 2 p'm' sec & 
& vm S cot a = C 8 T, i ■ >0 , 

a 2 ^ sec /3 — a' 2 p'm' sec /? 

cot(Xl <m 2 -Tm 3 sin 2 /? ' 

whence the change of phase by reflection is determined ; also, 

5 2 :5f:5' 2 
= (an 2 — Tm* sin 2 /3) 2 + (a\m sec /? + aPo'm' sec /3') 2 , 
:((r» 2 — Tm 2 sin 2 /?) + (a?om sec/3— a' 2 p'm' sec/?') 9 , 
: 4a 4 p 2 m 2 sec 2 /?', 

giving the ratios of the intensities of the incident, reflected and refracted waves. 
Put a = and T= and we obtain the results of the cases considered by 
Green in his paper on the Eeflection and Refraction of Sound, published in the 
Transactions of the Cambridge Philosophical Society, 1838, and republished by 
Ferrers in the Mathematical Papers of the late George Green, 1871. 
Then a = , a x = , 

and a P sec/9-aysec/y 

2ap 
D ap sec /9 -f- dp 1 sec /3' , 

jD = - O, 

2ap 
since am = a'ntf = n . 

Vol. IX. 



74 Gesenhill: Wave Motion in Hydrodynamics. 

We might, as in Green's paper, have supposed the incident, reflected and 
refracted plane waves given by 

<2>*= / { m (x sin /? — z cos /3) — nt + a \ , 
$ t = F { m (x sin /? + a cos /?) — r*# + a x j- , 
4>'=/i {m'(cc sin/3' — 3 cos/3')-— «£ + a'} , 
and the displacement of the curtain by 

n = / 2 |ma; sin /3 — nt) , 
when /, _F, /i , f % denote arbitrary functions, and determine the conditions to be 
satisfied as before. 

10. Theory of Long Waves in Canals. 

In this theory the vertical motion of the liquid particles is supposed insen- 
sible compared with the horizontal motion, and the depth of water small 
compared with the wave length ; so that 

U* = gh, 

as before in § 4, for water of uniform depth. 

This is proved independently by supposing the pressure at any depth the 
same as the hydrostatic pressure due to the depth below the free surface ; so that 
£ , denoting the horizontal displacement of a liquid particle, and r\ the elevation 

of the free surface, then 

d 2 f dp 

P !& ~ ~ ~dx 

.. dp dt) 

and *=W*i 

., . <?£ drj 

SOthat d?=-^&- 

But the equation of continuity leads to the condition 



6(A + <l + i) = M; 



or, v + h- = o, 

to one order of approximation, b denoting the breadth and h the depth of the 
canal of water. Then 

df - gh dx? > 
so that JJ—gh, 

U denoting the velocity of wave propagation. 
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We may, however, generalize this kind of motion, as was done by Kelland, 
by supposing the cross-section of the canal of any form, but uniform ; and then 
Kelland found IP = gA/b , 

where A is the area of the cross-section, and b, as before, the breadth at the surface 
of the water. 

For, in the more general case, where we consider the waves at the surface 
of separation of two liquids of densities o and o' filling a closed uniform hori- 
zontal pipe or conduit, so that A and A' denote the cross-sections of the pipe 
occupied by the liquids, and b the breadth of the plane of separation, we shall 



find 



TP = 



■JL 



JL 4- JL 

A ^ A 1 



at 



reducing to gA/b, where p' = 0. 

The simplest way to prove this is to suppose the motion made steady by 
applying the reversed velocity — U, equivalent to considering the motion rela- 
tive to an origin moving with velocity U. 

Then y, denoting the elevation of the surface, and u, u' the small additional 
velocities in the liquid due to the wave motion, 

{A + by,){U + v) = AU, 

(A> — br!)(U+u')z=A'U; 
or, Au + b Uvj = , 

A'vl — bUri=:0. 

Then, dp and dp', denoting the increments of pressure caused by the wave 
motion in the liquids just below and just above the surface of separation, 
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ay +gtf+ \ ? '(u+ wy- ±- 9 >u>=o; 

or, dp + gpyj+pUu = 0, 

dp'+ gp'y+p'Uu'=0. 

When there is no capillarity, etc., dp = dp', so that 

9 (p — p0»7 = (p'**'— H i7 

But, with a separating film of tension T, 

T% + &-& = *; 
so, if we assume that vi = a cos mx, 

we have T^\ = — Tm 2 )? 

aar 

= dp' — dp = g(p — p') >7 + (p«— p'«0 *7 

= <Kp-p'h-(-^ + ^)&^; 

When the liquids are bounded below and above by horizontal planes, at 
distances h and h! from the mean plane of separation, this equation becomes 



Tin* 



= (x + ¥) F! -^P-P'). 



an equation which will be found useful as a preliminary to the consideration of 
the Instability of Jets and its application to the flapping of sails and flags, investi- 
gated by Lord Rayleigh (Proceedings of the London Mathematical Society, Vol. 
X, No. 141). 

If the upper liquid had been moving with mean velocity V different to U, 
the preceding equations would be replaced by 

Au + b Tin = , 
A'u' — bU'n = o, 
dp + gpy + pUu = 0, 

dp' + gpri + p'U'u'z=:0, 

and then Tm*= (p i£ + p'^) b-g(p-p>). 
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11. Waves at the Surface of Separation of Two Liquids. 

The preceding case suggests the consideration of the general case of waves 
at the surface of separation of two liquids of different densities, and conse- 
quently a horizontal plane, when the liquids are either still except for the wave 
motion, or are flowing across each other with given mean uniform velocities, in 
which case the liquids must be bounded above and below by horizontal plane 
barriers if these velocities are not in the same direction. 
First, when the liquids are still, we must have 

<£> = A cosh m {z + h) cos {mx — nt) 
in the lower liquid, of depth h , as before, and 

4/ = A' cosh m {z — h') cos {mx — nt) 
in the upper liquid, of depth h', suppose ; and then, if 

yi = a sin {mx — nt) 
represents the displacement of the surface of separation, we must have 

d>l _ dfy _d$_ 
dt cte dz ' 
when z = 0, in order that there should be no separation of the liquids; conse- 
quently, — na = mA sinh mJi = — m A' sinh mh'. 
Again, from the hydrodynamical equation 

*.+,. + £ = * 

we obtain, at the surface of separation, 

9p + m + P -£■ = 

in the lower liquid, just below the surface of separation, and 

W + #1 + ^ = 
in the upper liquid, just above. 

Neglecting capillarity, etc., dp = dp', and, therefore, 

= ( — onA cosh mh + p'nA' cosh mh') cos {mx — nt) 

p — coth mh + p' — coth mh' ) n : 
r m r m J " 

so that ^ = 4 = A ' f~< , . 

mr m p coth mh -J- />' coth mh' 
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(Stokes, On the Theory of Oscillatory Waves, Cam. Phil. Trans., Vol. VIII, p. 
441 ; republished in Mathematical and Physical Papers, Vol. I, p. 212.) 
Putting p' = , we obtain 

XP = ^tanh mh = ^tanh ^, 
as before. 

When 1 is small compared with h and M, then mh and mh' are large, and 
we may replace coth mh and coth mh' by unity, and then 

7-7-2 _ 9_ P — P f 
^ *~~ ' i / 1 

which, when the ratio p'/p is small, as is the case of air on water, can be replaced 

b ^ *7>=^-ri-2^Y 

12. Suppose, now, the upper liquid is moving, like the wind, over the 
surface of the lower liquid with velocity V, and we wish to determine XT, the 
velocity of propagation of waves of length "k at the surface of separation. 

For generality, we shall suppose the lower liquid also moving with velocity 
V, and seek to determine the new relation connecting V, V and XT. 

The simplest way is to take a moving origin or plane of yz, moving with 
velocity V in the direction of the axis of x, the direction of wave propagation, 
and to consider the relative motion of the liquid, which will now be steady 
relatively to the moving co-ordinate axes. 

This is equivalent to supposing the motion made steady by impressing the 
reversed velocity — XT on the system. Then we must put 

4» = ( V — XT) x + A cosh m (z + h) cos mx, 
<£>' = ( V — XT) x + A' cosh m(z — h') cos mx ; 

so that, now introducing the conjugate current functions -^ and ij/, 
■>|>=(F — XT)z — A sinh m (z — h) sin mx, 
4 , ' == ( V — XT) z — A' sinh m (z — /*') sin mx . 
For the liquids not to separate, we must have, when z = , 4" = 4-' ; so that 

A sinh mh = — A' sinh mh' , 

the condition obtained otherwise by putting, whenz= 0, 
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Now, supposing that the displacement of the surface of separation is given 
by >7 = a sin mx , 

we must have, when z = , 

■^ = ( V — U)(t! — a sin mx) , 
4 1 ' = (V — U)(vi — a sin nix) ; 
so that A sinh mh= ( F — Z7) a, 

.A' sinh mh'= — ( V— TJ) a. 

Also, from the hydrodynamical equations, with the same notation as before, 

or, 3p + fl^ — pw J. ( V — TJ) cosh mA sin mx = , 

^P' + ^p'^ — p'»i J.' ( V — TJ) cosh mh! sin ma; = ; 

or, dp = — #p>? + pwi ( F — TJ) 2 coth wA >7 , 

dp' = — go'n — p'ni ( V -i- TTf coth mh! *i . 

If there is no capillarity, etc., at the surface of separation, dp = dp' ; so that 
g (p — p') = pm ( F— #) 2 coth mA + p'm ( V — ?7) 2 coth mtf, 

whence £7 is determined, when F, V, p, p' and A, A' are given. 

We have supposed here that the current velocities F and V are in the same 
direction as TJ, the wave velocity ; but if F and V make angles a and a! with TJ, 
then, in the above expression, F and F' must be replaced by F cos a and 
V cos a', the components F sin a and V sin a' of the currents perpen- 
dicular to the direction of propagation of the waves having no effect upon 
the determination of TJ (Encyclopaedia Britannica, 9th edition, article Hydro- 
mechanics). 

13. In the most general case, where the surface of separation is endowed 
with tension T, superficial density <r and flexural rigidity L, the condition to be 
satisfied at this surface is 

A*-n sffi-n sHy> 

dp'. 

Now, if 



eft? 


T*3- 

dx 2 


L d^ 


+ dp 


V = 


a sin mx, 




df ~~ 


vAz 

dx 2 


— — n 


V 


da? 


— m\, 


A _ 

dx'~ 


■ m\; 
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also, 

dp — dp'= — g{o — p')>7+ {pm(V— Uf cothmh + p'm(V — Uf coth mh' } 17 ; 

so that, omitting the common factor r\ , 

cm 4 — Tm 2 — Lm* — g(p — p') + pm( V— Uf coth mh + p'wi ( V— Uf coth m# = . 

The application of this equation to the discussion of the Instability of Jets, 
including the flapping of flags and sails, has been considered by Eayleigh, as 
mentioned above. 

In the application to a flag we may put p = p', and replace coth mh and 
coth mh' by unity ; also, L — 0, and we may also suppose T=- 0; then 

an 2 + 2pm (V— Uf=0, 

indicating the instability of the motion, and showing that it cannot be repre- 
sented by a periodic term of small displacement ; we must therefore replace in 
the motion cos (mx — nt) by cosh (mx — nt) , sinh (mx — nt), or 

(P cosh mx + Q sinh mx) cos nt. 

In the above general equation put <r = 0, i=0, F= 0, and replace 
coth mh and coth mh' by unity; then 

Tm* + g (p — p') = pmU 2 + p'm (V — Uf, 

the equation considered by Thomson for the determination of the ripples pro- 
duced by wind V over the surface of still water. 

If W is the velocity of ripples of the same wave length without wind, 
Tm 2 + g{p - p') = {p + p')mW 2 ; 

W 2 = ^^ l+ - mT 



m p + P ' ' p + {>" 
the minimum value of which, for different values of m , is 

W^2^gT)^=p, 
and then m* = g (p — p')/ T. 

But 9 U>+ f '(r— U)'=( 9 + f')W; 

so that U= IIL, ± /[ W , H>y\ 



V{ 



p + p' V { {p+yfl ' 

giving the velocities of the ripples with and against the wind V. 

The least value of V' 2 is less than (p + p')Vpp' times the least value of W 2 , 
and, therefore, the least value of V' 2 to produce ripples is 

W(gT) p -±£s/(p-p'). 

rr 
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If the wind is blowing with velocity greater than this minimum value of V, 
the plane surface of the water becomes unstable, and ripples are produced 
(Sir W. Thomson, Phil. Mag., 1871). 

14. Waves in Water Flowing with Variable Velocity fz, Some Function of the Depth z. 

In this manner we may attempt an investigation of the standing waves seen 
in a sloping current of water, where the velocity varies with the depth in conse- 
quence of viscosity and the fluid friction against the bottom ; the method, 
however, will not be very rigorous, as we must begin by assuming fluid friction 
to account for the varying velocities at different depths, and afterwards neglect 
fluid friction when we come to consider the superposed wave motion. 

Supposing, however, the motion is steady, we must put the current function 

4/ =. Fz — A sinh m (z + A) sin mx , 

so that the mean value of u or ~ is F'z or fz, denoting fz by F'z. Then, at 

the bottom of the water, 

4> = Fh , a constant, 

and at the surface 

^ = FO + vfO — A sinh mh sin mx ; 

so that if, at the surface, 

yj = a sin mx , 

then A sinh mh = afO . 

At the surface 

ai> + ^ + ip(§)-|p(/0) 2 = 0; 

or, dp + gori + p (jz/'O — mA cosh mh sin mx)fO = ; 

or, dp + or\ \g + /O/'O — m (fOf coth mh \ = ; 

so that, if dp = , 

9 + /O/'O — m (fOf coth mh = . 

Here /0 denotes the velocity of the current at the surface, and /'0 the vertical 
rate of change of the velocity at the surface. 

For instance, if the current flows uniformly with velocity V, 

g — m F 2 coth mh = ; 

or, T* = — tanhwiA, 

m 

as before. 

Vol. IX. 
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For a viscous liquid, flowing over a flat bottom, f"z = 0, f'z = 
denoting the viscosity and i the slope of the stream, supposed small, 



91 
¥■ 



V 






z + C, 



supposing V the current velocity at the surface (Bncy. Brit., Hydraulics). There- 
fore, g + VG — mV* coth mh = ; 

G is generally determined from the condition that the liquid adheres to the 
bottom, and, therefore, /& = 0, giving 

2 [i h 

15. In the experimental verification of the above theory of the motion of 
waves at the surface of separation of two liquids, we can make the wave velocity 
U as small as we please by making p and p' nearly equal. 

Again, in order to study experimentally the waves in water of uniform 
depth, the best plan to obtain uniformity of depth is to pour water on the top of 
mercury (Stokes, Math, and Phys. Papers, I, p. 199). But in this case the mer- 
cury forming the bottom of the water will not be fixed, but will itself be set 
into wave motion, and the modification thus introduced is considered by Stokes 
on p. 217. This is a particular case of the general conditions to be satisfied 
when waves are propagated at the surfaces of separation of a number of superin- 
cumbent liquids forming horizontal strata, and limited above and below by fixed 
horizontal planes. If the upper surface is free, the density of the highest stratum 
of liquid must be supposed zero. 

This general theorem has been worked out by Mr. R. R. Webb, and we 
shall proceed to investigate his results, which were given in the Math. Tripos 
Examination at Cambridge in Jan., 1884, as follows : 

A rectangular pipe whose faces are horizontal and vertical planes is com- 
pletely filled with n + 1 liquids ; show that the velocities v of propagation of 
waves of length % at the surfaces of separation of the strata are given by the 
equation 



= 0. 
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where A r = 27tv % /X (p r+ i- coth 2nh r+1 /h + p r coth 2nh r /X) — g (p r +i — p r ) , 

B r = 2nv % p r /% cosech 2nh r /X , 
and h r is the equilibrium thickness of the r th stratum of density p r . 

In particular, if p r — ra, h r = ra, then the In values of v are included in 
the formula ,l//\ / r n \ 

where r is supposed to assume the values 1, 2, 3 n, and %, the wave length, 

is supposed very large compared with na. 

In order that the velocity function ty r in the r th stratum should satisfy the 
equation of continuity and the conditions that 

9r _ 9>--i a t t]j e sur f ac e of separation of the r th and r — 1 th liquids, 

then, taking the plane of xy in the upper surface of the pipe, we must have 
<p r = I C r _ 1 cosh m (z + \ + h % + . . . . + h r ) 

- G r co S hm(z + h 1 + h % + - ■ ■ .+*,-!)} C -^^ . 
and then 

•>p = Cy_! cos (?na; — wtf) , when 2 = — h x — h 2 — .... — K-i, 

-Si = C r cos (rnx — nt) , when z = — h x — h % — . . . . — 7i r ; 

and, therefore, if v\ r denotes the elevation of the surface of separation between 
the r th and r + 1 th liquids, 

■— = G r cos (rnx — nt) , 

and ri r = sin (mx — nt) . 

To express the fact that there is no discontinuity at this surface of separa- 
tion, we have the equations 

£ + »* + & = •. 



i'rA-l 

so that, since dp r = dp r +\, 



— T y>7r -f 3^ 



3*r 1 3^) r+ i. 

gprVlr + Pr-df = OPr + l^r + pr+1 Q f ) 

.9 (Pr + X - Pr) ^ + 9 r + 1 ^ - p, ^ = 
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Now, when % = — \ — A 2 — . . . . — h r , 

-J£r = — ( <X_i cosech ni7i r — G r coth ni7i r ) sin (mx — ni) , 
at m x 

X +1 = — (G r coth mh r+1 — r+1 cosech mh r+1 ) sin (mx — nt) ; 

so that, dropping the common factor — sin (mx — nt) , 

— g(p r +i — p r ) G r -\ (Orp r +i coth mh r+1 — G r+1 p r+1 cosech m7i r+1 

— Cv-ipr cosech mh r + G r p r coth m7i,) = ; 
or, — G r _ x B r + G r A r — G r+1 B r+1 = 0, 

since n % /ra = 27tv % /%, ni = InfX. 

Also, the top and bottom of the pipe being fixed horizontal planes, G = 
and C„ + 1 = 0, so that the elimination of the C's leads to the determinant given 
above. 

When the pipe is open at the top, we can represent the motion by supposing 
pj = o , or B x = , and then the particular case for waves when there are two 
superincumbent liquids with a free surface has been given by Stokes (Math, 
and Phys. Papers, I, p. 217). 

It will be noticed that, although there is no discontinuity in the value of 

-¥- at a surface of separation of two strata, there is discontinuity in J±- , denoting 

a slipping of one surface over the other, the slipping velocity at the r th surface, 
where z = — h x — h % — .... — 7i r , being 

= { C r _ x cosech mh r — C,.(coth7nA r+1 +coth mh r )+ G r+1 cosech mh r + 1 \sm(mx — nt). 

This slipping proves a difficulty in the attempt of proceeding from the above 
investigation of waves in strata of finite thickness to the case of waves in a 
liquid of variable density arranged in horizontal strata. 

When h is large and in is consequently small, we may replace p r coth m7i r 
and p r cosech mh r by p r /mh r ; and then, if p r = ra, h r = ra, 

7l 2 

A r = 2cr —3 gaa = 2crw 2 — gaa , 

B r = a — r = av % : 
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so that the above determinant becomes 
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n rows 



ga 



= , where G— 2 — ^V , 

the determinant considered by Kayleigh (Theory of Sound, I, p. 131). 
It is there proved that, putting (7=2 cos $■, so that 



1 1 

v i = —gased i ^-^, 



then 
where 



rn 



r=l, 2, 3, 



n. 



1 6. Waves in Canals with Sloping Sides, or against a Sloping Beach, 

So far, the wave motion considered has only involved two co-ordinates, x 
and z, and might be considered limited by any two fixed vertical planes per- 
pendicular to the axis of y. 

In the case of the canal of uniform arbitrary cross-section, Kelland obtained 
the equation, JP = gA/b, 

for U, the wave velocity of long waves moving along the canal. 

Kelland, however, was successful in obtaining an exact expression for the 
motion of progressive waves in a straight canal the sides of which sloped down 
uniformly to an edge at an angle of 45°; he found that, taking the axis of x 
along this edge, we can put 

<?> = A cosh my cosh mz cos V 2 (mx — nt) , 
satisfying the equation of continuity 



da? 






and also the boundary condition ^ = ; 



or, 



c® dy d% 
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I, m, n denoting the direction cosines of the normal to the boundary. In this 
case the boundary conditions are 

^-^ = 0, when 3,-3 = 0, 

dy dz v 

Q + d$=0, when j, + « = <>, 

dy -dz * 

which are immediately seen to be satisfied. 
At the surface of the water, where z = h, 



I 






<p, 



for all values of x and y ; or, 

nil sinh mh = cosh mh, or ml = coth mh, 



and 



(J*=~=9=l- tanh mh, 
m, Itw m 



the same as for waves in water of depth h, but now 27t/X = ms/1. 

By transferring the axis of a; to the edge of the water on one bank, we 
obtain <p = A cosh m (y + h) cosh m (z — h) cos V 2 (mx — nt) , 

which, when h is made indefinitely great, can be replaced by 

<£ = Be-™^-* cos */2 (mx — nt), 
giving the motion for waves moving parallel to a shore sloping at 45°, the crests 
of the waves being perpendicular to the shore. 

This may easily be generalized, as Stokes has shown (Report on Recent 
Researches in Hydrodynamics, Math, and Phys. Papers, I, p. 167), for a shore 
sloping at any angle a, by putting 

4> = Be~ m { y <*>*«-* sin «) cos (mx — nt) , 




satisfying the equation of continuity and the boundary condition 

^2 sin a + ^5 cos a = , when y sin a + z cos a = ; 

dy dz 
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also, at the free surface % = , 



l~$- = <p; or, ml sin a = 1 ; 



so that Z7 2 = — = ^ = -£- sin a = ■£- sin a , 



n 2 g g . qX 

wr to" 1 m 2?r 

a denoting here the slope of the shore to the horizon. 

Analogous to Kelland's previous solution, we might put 
$ = A sinh my sinh mz sin V 2 (»na; — nt) , 

satisfying the equation of continuity and the boundary conditions, and at the free 

surface z = h, l-£ = d> gives 

ml = tanh /nA ; 

g 
or, C/ -2 = — coth w?A . 

m 

The shape of the free surface will be different in the two cases ; in the first, 

Kelland's case, when z = h, 

-~- = ^J- = m A sinh mA cosh my cos \/ 2 (mx — nt) ; 

so that ri = — sinh otA cosh wy sin V 2 («w; — n£) 

or the form v\ = a cosh my sin \/ 2 (ma; — w£) , 
an even function of y ; and in the second case, 

-J- = ->r- = mA cosh m7i sinh ot?/ sin \f 2 (»j;r — nt) ; 

so that the free surface is of the form 

57 = a sinh to?/ cos V 2 (mx — nt) , 
an odd function of y . 

Introducing capillarity on the free surface, but neglecting its effect at the 

contact of the surface with the bank, then the equation 

dz gp dz 3 

Tim 9 

gives ml sinh mh -\ sinh mh = cosh mh ; 

99 

or, ml H = coth mh ; 

99 

and, for the second kind of motion, 



ml H = tanh mh. 

99 
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17. Let us apply Kelland's expression to determine the progressive waves 
at the surface of separation of two liquids, each half filling a pipe, of which the 
cross-section is a square with a vertical diagonal, of length 2h . 

Taking diagonals of the square as axes of y and z , then we can put 
<£> = A cosh my cosh m (z + h) cos V 2 (mx — nt) , 
$>' = — A cosh my cosh m(z — h) cos V 2 (mx — nt) , 

satisfying the equation of continuity and the boundary conditions, except just 
where the surface of separation meets the boundary, the disturbing effect of 
which we shall neglect, although of course disturbing waves would be generated 
thereby. 

Then, as before, if 77 denotes the elevation of the surface of separation 
where the mean value of z = , 



— ^- = ^2 = mA cosh my sinh mh cos V 2 (mx 
eft 3s a v 



nt); 



so that 



mA 



also, 



y\ = -jx cosh my sinh mh sin \/ 2 (mx — nt) ; 

-i = wJV 2 cosh my cosh mh sin V 2 (mx — nt) = — 2 — n coth m7t , 



and 



eft 

#' _ 
eft 



= 2 — 57 coth m?A . 
m 
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Also, as before, 



or, dp + gori — 2p — 57 coth mA = , 

9p' + gp'y; + 2p' — y\ coth mA = ; 



so that if dp = dp', 



9 (p — p') = 2 — (p + P') coth m7*; 

M 2 ] ^ jO 

m s 2 m p -\- p' 



or, Z7 2 = -g = — -^ r , f , tanh m/t. 

Similarly, by putting 



d> = J. sinh jny sinh m (z + A) sin V 2 (»nx — nt) , 
$' = A sinh my sinh ?n (s — h) sin V 2 (mx — nt) , 
we should obtain 

U*=l.± P^Zf^ co th mh. 

4 m p -j- p 

18. Standing Waves across a Rectangular Channel. 

We shall find that, by replacing the hyperbolic functions of y and 2 partly 
by circular functions in the above solutions for progressive waves along the 
channel, we shall be able to solve the motion of standing waves in which the 
crests are parallel to the axis of the canal or channel. For, if we put 

d> = A (cos my cosh mz + cosh my cos mz) cos nt , 
or <p = A (sin my sinh mz + sinh my sin mz) sin nt, 

expressions which are equivalent to those obtained by Kirchoff (Ueber stehende 
Schwingungen einer schweren Fliissigkeit, Gesammelte Abhandlungen, II, p. 
440), we shall satisfy the equation of continuity and the boundary conditions 

Sy-&=0,wheny-z=0, 

^ + d ± = 0,wheny + z = 0; 

dy dz 

and at the free surface z = h the condition 

'§? = *- 



90 Gebenhill : Wave Motion in Hydrodynamics. 

for all values of y, leads to the equations 

ml (cos my sinh mh — cosh my sin mh) = cos my cosh mh + cosh my cos mh, 

giving ml sinh mh = cosh mh and — ml sin mh = cos ?«A, 
or ml = coth mA = — cot mh , 

for the even vibrations, and 

ml (sin my cosh «iA + sinh my cos wA) = sin my sinh wiA + sinh my sin mA , 

giving ml cosh mh = sinh mh and mZ cos mh = sin mA, 

or ml = tanh wiA = tan mh , 

for the odd vibrations. 

In these equations mh is the same as Kirchoff 's p, and, with the notation 

of the hyperbolic functions, KirchofFs period equations 

. , sin p 

sinh»= — — ±~- 
* V( cos 2p) 

__ cosp 



and cosh «= „ 

* V( cos *p) 

correspond to our period equations 

coth p = — cot jo 
and tanh j> = tan j» ; 

both being included in the single equation 

cos 2p cosh 2p = 1 , 

the period equation for the lateral vibrations of a free-free or clamped-clamped 
bar (Kayleigh, Sound, I, p. 219). 

In fact, the vibrations of the surface are of exactly the same character as 
those of a free-free bar of length 2h, the first value of q> giving the even and 
the second the odd vibrations. 

Suppose the surface were covered with ice of thickness e and flexural rigidity 

L = yr- $E) then, at the surface z = h, 



12 



^1 t diy } i o> 



and dp + gov; + p S = 0. 

But §=-n»n,^ = m* V ; 

go that rfperi = m i Lyj + got] + P -^f , 

(jt 
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or \p(g — n % e)+m'L}vi= — p£; 

and since * = IMS? = ~ "**■ 



or since g = In*, 



\o{g — n*e) + m i L\^ = n\$ ] 







and, therefore, as before, 

mil — e) -\ — 5— = coth rnh = — cot mh, 

or =r tanh rnh = tan mh , 

showing that the length of the equivalent simple pendulum is altered by 
m^.L/tfo — e by the presence of the ice. 

Dropping for the present the factors A cos nt and A sin nt, then 
d> = cos my cosh mz + cosh my cos mz , 
or $ = sin my sinh mz + sinh my sin mz , 

and, therefore, the conjugate current functions are 

4* = sin my sinh mz — sinh my sin mz , 
or 4 = — cos my cosh «iz + cosh my cos mz ; 

so that ^> + *4 = cos w? (z + %) + cosh m (z + *y) , 

or = £ cos m (z + %) — i cosh m (z + iy) . 

Denoting $ + ii^ by w, and z + iy by ^«, then 

<w> = cos mu + cos imu , 
or w = i cos «m — t cos »», 

gives the required motion in a rectangular channel. 

By transferring the axis of x to the edge of the water on one bank, we 
obtain $ ~ cos m {y — ^) cosh m (z + h) + cosh m (y — h) cos m (z + h) , 
or d> = sin m (y — h) sinh m (z + 7&) + sinh m(y — h) sin m (z + A) ; 

and these are the co-ordinates employed by Kirchoff. 

When h is made indefinitely great, these expressions may be replaced by 

(omitting constant factors, and remembering that the axis of z is drawn vertically 

upwards) , 

$ = e m * (cos my — sin my) + e~ my (cos mz — sin mz) , 

and then, when z = , 7 9<2> 

dz 
if ???Z = 1 , 
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giving Kirchoff 's solution of standing waves parallel to a shore sloping at 45° 
(Abhandlungen, II> p. 434). 

Returning to the original axes, the second value of $ gives, when m is small, 
so that w? may be neglected, 

q> = 2myz ; 

or, restoring the periodic factor, we may put 

4> = 2myz sin nt, 
and then l = h; 
also, dtj 3<£ 

so that 



dt =-£ = tony Bin nt; 



~ m 
n = — 2 — y cos nt, 




so that the free surface of the liquid remains plane during this kind of wave 
motion. Also, ^ — ni(y 2 — z*) sin nt, 

so that the liquid particles oscillate in rectangular hyperbolas (Kirchoff, Abhand- 
lungen, II, p. 436). 



19. Standing Waves across a Channel of 120°. 

Let us now attempt the solution of the corresponding waves in a canal the 
sides of which slope at angles of 30° to the horizon, and are therefore inclined 
to each other at an angle of 1 20°. 
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First, we notice that we can begin with an algebraical solution by putting, 

as the simplest case, 

<£= J.<J> cos nt, 

where 3> = z 3 — 3y*z + h 3 , 

and the corresponding stream function 

f = Zyi— y 3 
= y(zV3 — y)(zV3 + y), 

which vanishes when y = ± z\/3, showing that the boundary conditions are 

satisfied, and also, 

4> + M = (z + %) 3 + A 8 . 

Then, at the free surface z = h, 

<® = 3¥— 3f, <D = 2h 3 - 3% 3 ; 



so that 






and, thei'efore, l = h. 

The free surface is now a parabolic cylinder, for 

^- = ZA{h % ~ y*) cos nt, 

Civ 

yc=- — (W — iF) sin nt. 




For waves of a higher order, let us attempt the solution by putting 
w = cos tnu + cos rn(3u -j- cos m(!Pu, 
where m = 2 + iy, 

w = <& + i x 9, 

and /3 3 = — 1, $=\-\wZ, P~-\ \WZ. 
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Then, at the boundary y ■=.%*/%, 
w = cos (1 + £\/3) mz + cos -~- (1 — »V3) (1 + W3) mz + cos -5- (1 + i\Z3fmz 

= cos (1 + i\/3) mz + cos 2mz + cos ( — 1 + is/Z) z, 

a real quantity, so that I s = . 

Again, at the boundary y = — zV3 , w is real and ^ = 0. These conditions 
will also be satisfied by putting 

w = sin mu — sin m(3u + sin m^u , 

so that generally we can put 

w = sin m (u — a) — sin m ((3u -f- a) + sin m (fi^u — a) ; 
so that, since 

(3u + a = y (^V3 + z + 2a) + -y *' (y — ss*/3), 
£*w — a = — (y^/3 — z — 2a) 2" * (^ + S/ v /3 )> 

«i> = sin m (a — a) cosh wiy — sin -^ m (yVS + z -f- 2a) cosh-^- m (?/ — z\/3) 

+ sin -g- to (2/^3 — z — 2a) cosh -«r- to (y + z«/3) , 

ip = cos m (s — a) sinh my — cos -z- to (?/ V 3 + z + 2a) sinh -^ m(y — z V 3) 

— cosym(2/V3 — z — 2a) sinh -~ to (y + z\/3). 

Putting ?/ = z\/3, 

tp = cos to (z — a) sinh mz«/3 — cos m (z — a) sinh toz«/3 = , 

and putting y = — z\/ 3 , 

f = — cos m (z — a) sinh mz\/3 — cos m ( — z + a) sinh ( — mz\/3) = , 

so that the boundary conditions are satisfied. 

Expanded in ascending powers of (z •+- ix) , we shall find 

w = — 3 sin a — -5- m 3 (z + *'a5) 3 cos a + 

so that when m is small, and to 5 , ... . can be neglected, we obtain the previous 
algebraical solution. 

At the free surface z = h we must have 

as 
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for all values of y ; and, therefore, since we may write 

<3> = sin m (z — a) cosh my 

1 111 

— 2 sin -jr- m (z + 2a) cosh -y mzis/ 3 cosh -„- to?/ cos -~- rny\/ 3 

+ 2 cos -Q- m (z + 2a) sinh -~- mz */ 3 sinh -*- my sin -x- my\/ 3 , 

l~ ■=. ml cos to (z — a) cosh my 
dz 

— ml \ cos -x- to (s + 2a) cosh -x- mz\/ 3 + \/ 3 sin -_- to (z + 2a) 

sinh -5- mz/y/3 [ cosh -x- my cos -x- my^Z 

( 1 1 1 

+ mn — sin -x- to (z + 2a) sinh -_- toz\/ 3 + V 3 cos -x- m (z + 2a) 

cosh-x- mzs/2> \ sinh -x- to?/ sin -x- wz*/3 ; 

therefore, at the free surface z = h, 

ml cos m (h — a) = sin m (Ji — a) , (I) 

ml J cos -x- m (h + 2a) cosh -x- mh\/3 + \/3 sin -x- to (A -f 2a) sinh -x- mh*/3 I 

= 2 sin -~- to (A H- 2a) cosh -_- mh\/3 , (II) 

to?| — sin -x- to (A + 2a) sinh -_- mhV3+\/3 cos -x- to (7*+ 2a) cosh-x- mh\/3 [ 

= 2 cos yffl(H 2a) sinh -_- TO7i\/ 3 , (III) 

three equations for determining mh, a and I. 
From (II) and (III), eliminating ml, 

cos 2 -x- to (A + 2a) sinh mhs/Z + */3 sin to (A + 2a) sinh 2 -x- mh\/3 

= — sin 2 -x- to (7;. + 2a) sinh ?rrA V 3 -+- V 3 sin to (A + 2a) cosh 2 -x- toA V 3 ; 

or, sinh mh*/S = */3 sin to (A + 2a). (IV) 

Also, from (II) and (III), 

1 to? cosh -x- to7m/ 3 

tan -x to (A + 2a) = 



2 cosh -x- mhis/ 3 — to?>/ 3 sinh -x- w?7i\/ 3 

1 . 1 

ml^/Z cosh -x- mh/s/3 — 2 sinh -x- m7i>/3 

. - ; 

ml sinh -x- to7«\/3 
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so that 

■y m?P sinh mh^/3 = 2m?v'3 f eosh 2 -^ mA\/ 3 + sinh 2 --- mhs/Z) 

3 
;r- m 2 Z 2 sinh mhs/3 — 2 sinh mA\/3 ; 

or, m? + — ^ =\/3 coth mA\/3. (V) 

From (I), ml = tan m (ft — a) , 

so that 

ml + . — j = 2 cosec 2m (A — a) ; 

2 
or, sin 2m (A — a) = ~-jk tanh mh */?>. (VI) 

From (IV) and (VI), 

sin mA cos 2ma + cos mA sin 2ma = —^ sinh mh */3, 

2 
sin 2mA cos 2ma — cos 2mh sin 2ma = ~r? tanh mh \/3 ; 

and, therefore, 

<\/3 sin 3m7t cos 2ma = sinh mh^/3 (cos 2mA A = — ; - ) , 

\ cosh mh a/3/ 

\/3 sin 3mA sin 2ma = sinh mh V3 ( sin 2mA — — : — = ) . 

\ cosh mh \/3/ 

Squaring and adding, to eliminate a, 

3 sin 8 3 mh 4 cos 3 mh 4 

1 i il — l — 7^ + 



sinh' mh v^3 cosh mh \/3 cosh 8 mh \/3 ' 

or, denoting cos 3mA by a, and cosh mh \/3 by /?, 

* p — 1 ^ /3 T /? 8 ' 

reducing to 

/3 4 + 4a/3 3 + 3a 2 /3 2 — 4a/? — 4 = 0, 
or (/3 2 +2a/3) 2 = (a/3+2) 2 , 

/3 2 +2a/3 = ±(a/? + 2); 
so that /? 2 + a/3=2, 

or /3 2 + 3a/? = — 2 ; 

a=— + 2/73, 
or 3a =— /?— 2/0; 

cos 3mA = — cosh mh*/Z + 2 sech mA <\/3 , 
or 3 cos 3mA = — coshmAV3 — 2sechmAy'3, 

the period equations. 



